Multimodal anatomical mapping of subcortical regions in marmoset monkeys using high-resolution MRI and matched histology with multiple stains.

Subcortical nuclei and other deep brain structures play essential roles in regulating the central and peripheral nervous systems. However, many of these nuclei and their subregions are challenging to identify and delineate in conventional MRI due to their small size, hidden location, and often subtle contrasts compared to neighboring regions. To address these limitations, we scanned the whole brain of the marmoset monkeys in ex vivo using a clinically feasible diffusion MRI method, called the mean apparent propagator (MAP)-MRI, along with T2W and MTR (T1-like contrast) images acquired at 7 Tesla. Additionally, we registered these multimodal MRI volumes to the high-resolution images of matched whole-brain histology sections with seven different stains obtained from the same brain specimens. At high spatial resolution, the microstructural parameters and fiber orientation distribution functions derived with MAP-MRI can distinguish the subregions of many subcortical and deep brain structures, including fiber tracts of different sizes and orientations. The good correlation with multiple but distinct histological stains from the same brain serves as a thorough validation of the structures identified with MAP-MRI and other MRI parameters. Moreover, the anatomical details of deep brain structures found in the volumes of MAP-MRI parameters are not visible in conventional T1W or T2W images. The high-resolution mapping using novel MRI contrasts, combined and correlated with histology, can elucidate structures that were previously invisible radiologically. Thus, this multimodal approach offers a roadmap toward identifying salient brain areas in vivo in future neuroradiological studies. It also provides a useful anatomical standard reference for the region definition of subcortical targets and the generation of a 3D digital template atlas for the marmoset brain research (Saleem et al., 2023). Additionally, we conducted a cross-species comparison between marmoset and macaque monkeys using results from our previous studies (Saleem et al., 2021). We found that the two species had distinct patterns of iron distribution in subregions of the basal ganglia, red nucleus, and deep cerebellar nuclei, confirmed with T2W MRI and histology.

Multimodal anatomical mapping of subcortical regions in Marmoset monkeys using high-resolution MRI and matched histology with multiple stains.

Subcortical nuclei and other deep brain structures play essential roles in regulating the central and peripheral nervous systems. However, many of these nuclei and their subregions are challenging to identify and delineate in conventional MRI due to their small size, hidden location, and often subtle contrasts compared to neighboring regions. To address these limitations, we scanned the whole brain of the marmoset monkeys in ex vivo using a clinically feasible diffusion MRI method, called the mean apparent propagator (MAP)-MRI, along with T2W and MTR (T1-like contrast) images acquired at 7 Tesla. Additionally, we registered these multimodal MRI volumes to the high-resolution images of matched whole-brain histology sections with seven different stains obtained from the same brain specimens. At high spatial resolution, the microstructural parameters and fiber orientation distribution functions derived with MAP-MRI can distinguish the subregions of many subcortical and deep brain structures, including fiber tracts of different sizes and orientations. The good correlation with multiple but distinct histological stains from the same brain serves as a thorough validation of the structures identified with MAP-MRI and other MRI parameters. Moreover, the anatomical details of deep brain structures found in the volumes of MAP-MRI parameters are not visible in conventional T1W or T2W images. The high-resolution mapping using novel MRI contrasts, combined and correlated with histology, can elucidate structures that were previously invisible radiologically. Thus, this multimodal approach offers a roadmap toward identifying salient brain areas in vivo in future neuroradiological studies. It also provides a useful anatomical standard reference for the region definition of subcortical targets and the generation of a 3D digital template atlas for the marmoset brain research (Saleem et al., 2023). Additionally, we conducted a cross-species comparison between marmoset and macaque monkeys using results from our previous studies (Saleem et al., 2021). We found that the two species had distinct patterns of iron distribution in subregions of the basal ganglia, red nucleus, and deep cerebellar nuclei, confirmed with T2W MRI and histology.

An integrated resource for functional and structural connectivity of the marmoset brain.

Comprehensive integration of structural and functional connectivity data is required to model brain functions accurately. While resources for studying the structural connectivity of non-human primate brains already exist, their integration with functional connectivity data has remained unavailable. Here we present a comprehensive resource that integrates the most extensive awake marmoset resting-state fMRI data available to date (39 marmoset monkeys, 710 runs, 12117 mins) with previously published cellular-level neuronal tracing data (52 marmoset monkeys, 143 injections) and multi-resolution diffusion MRI datasets. The combination of these data allowed us to (1) map the fine-detailed functional brain networks and cortical parcellations, (2) develop a deep-learning-based parcellation generator that preserves the topographical organization of functional connectivity and reflects individual variabilities, and (3) investigate the structural basis underlying functional connectivity by computational modeling. This resource will enable modeling structure-function relationships and facilitate future comparative and translational studies of primate brains.

An open access resource for functional brain connectivity from fully awake marmosets.

The common marmoset (Callithrix jacchus) is quickly gaining traction as a premier neuroscientific model. However, considerable progress is still needed in understanding the functional and structural organization of the marmoset brain to rival that documented in longstanding preclinical model species, like mice, rats, and Old World primates. To accelerate such progress, we present the Marmoset Functional Brain Connectivity Resource (marmosetbrainconnectome.org), currently consisting of over 70 h of resting-state fMRI (RS-fMRI) data acquired at 500 µm isotropic resolution from 31 fully awake marmosets in a common stereotactic space. Three-dimensional functional connectivity (FC) maps for every cortical and subcortical gray matter voxel are stored online. Users can instantaneously view, manipulate, and download any whole-brain functional connectivity (FC) topology (at the subject- or group-level) along with the raw datasets and preprocessing code. Importantly, researchers can use this resource to test hypotheses about FC directly - with no additional analyses required - yielding whole-brain correlations for any gray matter voxel on demand. We demonstrate the resource's utility for presurgical planning and comparison with tracer-based neuronal connectivity as proof of concept. Complementing existing structural connectivity resources for the marmoset brain, the Marmoset Functional Brain Connectivity Resource affords users the distinct advantage of exploring the connectivity of any voxel in the marmoset brain, not limited to injection sites nor constrained by regional atlases. With the entire raw database (RS-fMRI and structural images) and preprocessing code openly available for download and use, we expect this resource to be broadly valuable to test novel hypotheses about the functional organization of the marmoset brain.

High-resolution mapping and digital atlas of subcortical regions in the macaque monkey based on matched MAP-MRI and histology.

Subcortical nuclei and other deep brain structures are known to play an important role in the regulation of the central and peripheral nervous systems. It can be difficult to identify and delineate many of these nuclei and their finer subdivisions in conventional MRI due to their small size, buried location, and often subtle contrast compared to neighboring tissue. To address this problem, we applied a multi-modal approach in ex vivo non-human primate (NHP) brain that includes high-resolution mean apparent propagator (MAP)-MRI and five different histological stains imaged with high-resolution microscopy in the brain of the same subject. By registering these high-dimensional MRI data to high-resolution histology data, we can map the location, boundaries, subdivisions, and micro-architectural features of subcortical gray matter regions in the macaque monkey brain. At high spatial resolution, diffusion MRI in general, and MAP-MRI in particular, can distinguish a large number of deep brain structures, including the larger and smaller white matter fiber tracts as well as architectonic features within various nuclei. Correlation with histology from the same brain enables a thorough validation of the structures identified with MAP-MRI. Moreover, anatomical details that are evident in images of MAP-MRI parameters are not visible in conventional T1-weighted images. We also derived subcortical template "SC21" from segmented MRI slices in three-dimensions and registered this volume to a previously published anatomical template with cortical parcellation (Reveley et al., 2017; Saleem and Logothetis, 2012), thereby integrating the 3D segmentation of both cortical and subcortical regions into the same volume. This newly updated three-dimensional D99 digital brain atlas (V2.0) is intended for use as a reference standard for macaque neuroanatomical, functional, and connectional imaging studies, involving both cortical and subcortical targets. The SC21 and D99 digital templates are available as volumes and surfaces in standard NIFTI and GIFTI formats.

Marmoset Brain Mapping V3: Population multi-modal standard volumetric and surface-based templates.

The standard anatomical brain template provides a common space and coordinate system for visualizing and analyzing neuroimaging data from large cohorts of subjects. Previous templates and atlases for the common marmoset brain were either based on data from a single individual or lacked essential functionalities for neuroimaging analysis. Here, we present new population-based in-vivo standard templates and tools derived from multi-modal data of 27 marmosets, including multiple types of T1w and T2w contrast images, DTI contrasts, and large field-of-view MRI and CT images. We performed multi-atlas labeling of anatomical structures on the new templates and constructed highly accurate tissue-type segmentation maps to facilitate volumetric studies. We built fully featured brain surfaces and cortical flat maps to facilitate 3D visualization and surface-based analyses, which are compatible with most surface analyzing tools, including FreeSurfer, AFNI/SUMA, and the Connectome Workbench. Analysis of the MRI and CT datasets revealed significant variations in brain shapes, sizes, and regional volumes of brain structures, highlighting substantial individual variabilities in the marmoset population. Thus, our population-based template and associated tools provide a versatile analysis platform and standard coordinate system for a wide range of MRI and connectome studies of common marmosets. These new template tools comprise version 3 of our Marmoset Brain Mapping Project and are publicly available via marmosetbrainmapping.org/v3.html.

A resource for the detailed 3D mapping of white matter pathways in the marmoset brain.

While the fundamental importance of the white matter in supporting neuronal communication is well known, existing publications of primate brains do not feature a detailed description of its complex anatomy. The main barrier to achieving this is that existing primate neuroimaging data have insufficient spatial resolution to resolve white matter pathways fully. Here we present a resource that allows detailed descriptions of white matter structures and trajectories of fiber pathways in the marmoset brain. The resource includes: (1) the highest-resolution diffusion-weighted MRI data available to date, which reveal white matter features not previously described; (2) a comprehensive three-dimensional white matter atlas depicting fiber pathways that were either omitted or misidentified in previous atlases; and (3) comprehensive fiber pathway maps of cortical connections combining diffusion-weighted MRI tractography and neuronal tracing data. The resource, which can be downloaded from marmosetbrainmapping.org, will facilitate studies of brain connectivity and the development of tractography algorithms in the primate brain.

Anatomical and functional investigation of the marmoset default mode network.

The default mode network (DMN) is associated with a wide range of brain functions. In humans, the DMN is marked by strong functional connectivity among three core regions: medial prefrontal cortex (mPFC), posterior parietal cortex (PPC), and the medial parietal and posterior cingulate cortex (PCC). Neuroimaging studies have shown that the DMN also exists in non-human primates, suggesting that it may be a conserved feature of the primate brain. Here, we found that, in common marmosets, the dorsolateral prefrontal cortex (dlPFC; peak at A8aD) has robust fMRI functional connectivity and reciprocal anatomical connections with the posterior DMN core regions (PPC and PCC), while the mPFC has weak connections with the posterior DMN core regions. This strong dlPFC but weak mPFC connectivity in marmoset differs markedly from the stereotypical DMN in humans. The mPFC may be involved in brain functions that are further developed in humans than in other primates.

Neuroprotective Effects of MAGL (Monoacylglycerol Lipase) Inhibitors in Experimental Ischemic Stroke.

BACKGROUND AND PURPOSE: MAGL (monoacylglycerol lipase) is an enzyme that hydrolyzes the endocannabinoid 2-arachidonoylglycerol and regulates the production of arachidonic acid and prostaglandins-substances that mediate tissue inflammatory response. Here, we have studied the effects of the selective MAGL inhibitors JZL184 and MJN110 and their underlying molecular mechanisms on 3 different experimental models of focal cerebral ischemia. METHODS: SHR (spontaneously hypertensive rats) and normotensive WKY (Wistar Kyoto) rats were subject to an intracortical injection of the potent vasoconstrictor endothelin-1, permanent occlusion of a distal segment of the middle cerebral artery via craniectomy, or transient occlusion of the middle cerebral artery by the intraluminal suture method. JZL184 or MJN110 was administered 60 minutes after focal cerebral ischemia. Infarct volumes, hemispheric swelling, and functional outcomes were assessed between days 1 to 28 by magnetic resonance imaging, histology, and behavioral tests. RESULTS: Pharmacological inhibition of MAGL significantly attenuated infarct volume and hemispheric swelling. MAGL inhibition also ameliorated sensorimotor deficits, suppressed inflammatory response, and decreased the number of degenerating neurons. These beneficial effects of MAGL inhibition were not fully abrogated by selective antagonists of cannabinoid receptors, indicating that the anti-inflammatory effects are caused by inhibition of eicosanoid production rather than by activation of cannabinoid receptors. CONCLUSIONS: Our results suggest that MAGL may contribute to the pathophysiology of focal cerebral ischemia and is thus a promising therapeutic target for the treatment of ischemic stroke.

Investigating the spatiotemporal characteristics of the deoxyhemoglobin-related and deoxyhemoglobin-unrelated functional hemodynamic response across cortical layers in awake marmosets.

Blood oxygenation level dependent (BOLD) functional magnetic resonance imaging (fMRI) has become a major tool to map neural activity. However, the spatiotemporal characteristics of the BOLD functional hemodynamic response across the cortical layers remain poorly understood. While human fMRI studies suffer from low spatiotemporal resolution, the use of anesthesia in animal models introduces confounding factors. Additionally, inflow contributions to the fMRI signal become non-negligible when short repetition times (TRs) are used. In the present work, we mapped the BOLD fMRI response to somatosensory stimulation in awake marmosets. To address the above technical concerns, we used a dual-echo gradient-recalled echo planar imaging (GR-EPI) sequence to separate the deoxyhemoglobin-related response (absolute T2* differences) from the deoxyhemoglobin-unrelated response (relative S0 changes). We employed a spatial saturation pulse to saturate incoming arterial spins and reduce inflow effects. Functional GR-EPI images were obtained from a single coronal slice with two different echo times (13.5 and 40.5ms) and TR=0.2s. BOLD, T2*, and S0 images were calculated and their functional responses were detected in both hemispheres of primary somatosensory cortex, from which five laminar regions (L1+2, L3, L4, L5, and L6) were derived. The spatiotemporal distribution of the BOLD response across the cortical layers was heterogeneous, with the middle layers having the highest BOLD amplitudes and shortest onset times. ΔT2* also showed a similar trend. However, functional S0 changes were detected only in L1+2, with a fast onset time. Because inflow effects were minimized, the source of S0 functional changes in L1+2 could be attributed to a reduction of cerebrospinal fluid volume fraction due to the functional increase in cerebral blood volume and to unmodeled T2* changes in the extra- and intra-venous compartments. Caution should be exercised when interpreting laminar BOLD fMRI changes in superficial layers as surrogates of underlying neural activity.

A digital 3D atlas of the marmoset brain based on multi-modal MRI.

The common marmoset (Callithrix jacchus) is a New-World monkey of growing interest in neuroscience. Magnetic resonance imaging (MRI) is an essential tool to unveil the anatomical and functional organization of the marmoset brain. To facilitate identification of regions of interest, it is desirable to register MR images to an atlas of the brain. However, currently available atlases of the marmoset brain are mainly based on 2D histological data, which are difficult to apply to 3D imaging techniques. Here, we constructed a 3D digital atlas based on high-resolution ex-vivo MRI images, including magnetization transfer ratio (a T1-like contrast), T2w images, and multi-shell diffusion MRI. Based on the multi-modal MRI images, we manually delineated 54 cortical areas and 16 subcortical regions on one hemisphere of the brain (the core version). The 54 cortical areas were merged into 13 larger cortical regions according to their locations to yield a coarse version of the atlas, and also parcellated into 106 sub-regions using a connectivity-based parcellation method to produce a refined atlas. Finally, we compared the new atlas set with existing histology atlases and demonstrated its applications in connectome studies, and in resting state and stimulus-based fMRI. The atlas set has been integrated into the widely-distributed neuroimaging data analysis software AFNI and SUMA, providing a readily usable multi-modal template space with multi-level anatomical labels (including labels from the Paxinos atlas) that can facilitate various neuroimaging studies of marmosets.

Design and implementation of embedded 8-channel receive-only arrays for whole-brain MRI and fMRI of conscious awake marmosets.

PURPOSE: The common marmoset (Callithrix jacchus) is a New World primate of increasing interest to neuroscience and in translational brain research. The present work describes the design and implementation of individualized 8-channel receive-only radiofrequency (RF) coil arrays that provide whole-brain coverage and allow anatomical and functional MRI experiments in conscious, awake marmosets. METHODS: The coil arrays were designed with their elements embedded inside individualized restraint helmets. The size, geometry, and arrangement of the coil elements were optimized to allow whole-brain coverage. Coil-to-coil decoupling was achieved by a combination of geometric decoupling and low input impedance preamplifiers. The performance of the embedded arrays was compared against that of one 8-channel receive-only array built to fit the external surface of the helmets. RESULTS: Three individualized helmets with embedded coil arrays were built for three marmosets. Whole-brain coverage was achieved with high sensitivity extending over the entire cortex. Visual stimulation of conscious awake marmosets elicited robust BOLD fMRI responses in both primary and higher order visual areas of the occipitotemporal cortex. CONCLUSION: The high sensitivity provided by embedded receive-only coil arrays allows both anatomical and functional MRI data to be obtained with high spatial resolution in conscious, awake marmosets. Magn Reson Med 78:387-398, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Functional Connectivity Hubs and Networks in the Awake Marmoset Brain.

In combination with advances in analytical methods, resting-state fMRI is allowing unprecedented access to a better understanding of the network organization of the brain. Increasing evidence suggests that this architecture may incorporate highly functionally connected nodes, or "hubs", and we have recently proposed local functional connectivity density (lFCD) mapping to identify highly-connected nodes in the human brain. Here, we imaged awake nonhuman primates to test whether, like the human brain, the marmoset brain contains FC hubs. Ten adult common marmosets (Callithrix jacchus) were acclimated to mild, comfortable restraint using individualized helmets. Following restraint training, resting BOLD data were acquired during eight consecutive 10 min scans for each subject. lFCD revealed prominent cortical and subcortical hubs of connectivity across the marmoset brain; specifically, in primary and secondary visual cortices (V1/V2), higher-order visual association areas (A19M/V6[DM]), posterior parietal and posterior cingulate areas (PGM and A23b/A31), thalamus, dorsal and ventral striatal areas (caudate, putamen, lateral septal nucleus, and anterior cingulate cortex (A24a). lFCD hubs were highly connected to widespread areas of the brain, and further revealed significant network-network interactions. These data provide a baseline platform for future investigations in a nonhuman primate model of the brain's network topology.

An embedded four-channel receive-only RF coil array for fMRI experiments of the somatosensory pathway in conscious awake marmosets.

fMRI has established itself as the main research tool in neuroscience and brain cognitive research. The common marmoset (Callithrix jacchus) is a non-human primate model of increasing interest in biomedical research. However, commercial MRI coils for marmosets are not generally available. The present work describes the design and construction of a four-channel receive-only surface RF coil array with excellent signal-to-noise ratio (SNR) specifically optimized for fMRI experiments in awake marmosets in response to somatosensory stimulation. The array was designed as part of a helmet-based head restraint system used to prevent motion during the scans. High SNR was obtained by building the coil array using a thin and flexible substrate glued to the inner surface of the restraint helmet, so as to minimize the distance between the array elements and the somatosensory cortex. Decoupling between coil elements was achieved by partial geometrical overlapping and by connecting them to home-built low-input-impedance preamplifiers. In vivo images show excellent coverage of the brain cortical surface with high sensitivity near the somatosensory cortex. Embedding the coil elements within the restraint helmet allowed fMRI data in response to somatosensory stimulation to be collected with high sensitivity and reproducibility in conscious, awake marmosets.

BOLD responses to different temporal frequency stimuli in the lateral geniculate nucleus and visual cortex: insights into the neural basis of fMRI.

The neural basis of the blood oxygenation level dependent (BOLD) functional magnetic resonance imaging (fMRI) remains largely unknown after decades of research. To investigate this issue, the unique property of the temporal frequency tuning that could separate neural input and output in the primary visual cortex was used as a model. During moving grating stimuli of 1, 2, 10 and 20Hz temporal frequencies, we measured 9.4-T BOLD fMRI responses simultaneously in the primary visual cortex of area 17 (A17) and area 18 (A18), and the lateral geniculate nucleus (LGN) of isoflurane-anesthetized cat. Our results showed that preferred temporal frequencies of the BOLD responses for A17, A18 and LGN were 3.1Hz, 4.5Hz and 6.0Hz, respectively, which were comparable to the previously reported electrophysiological data. Additionally, the difference of BOLD response onset time between LGN and A17 was 0.5s, which is 18 times larger than the difference of neural activity onset time between these areas. We then compared the frequency-dependent BOLD fMRI response of A17 with tissue partial pressure of oxygen (pO(2)) and electrophysiological data of the same animal model reported by Viswanathan and Freeman (Nature Neuroscience, 2007). The BOLD tuning curve resembled the low frequency band (<12Hz) of local field potential (LFP) tuning curve rather than spiking activity, gamma band (25-90Hz) of LFP, and tissue pO(2) tuning curves, suggesting that the BOLD fMRI signal relates closer to low frequency LFP.